Order and Disorder in AKLT Antiferromagnets in Three Dimensions 
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The models constructed by Affleck, Kennedy, Lieb, and Tasakji describe a family of quantum 
antiferromagnets on arbitrary lattices, where the local spin S is an integer multiple M of half the 
lattice coordination number. The equal time quantum correlations in their ground states may be 
computed as finite temperature correlations of a classical 0(3) model on the same lattice, where 
the temperature is given by T = 1/M. In dimensions d = 1 and d — 2 this mapping implies that 
all AKLT states are quantum disordered. We consider AKLT states in d = 3 where the nature 
of the AKLT states is now a question of detail depending upon the choice of lattice and spin; for 
sufflciently large S some form of Neel order is almost inevitable. On the unfrustrated cubic lattice, 
we find that all AKLT states are ordered while for the unfrustrated diamond lattice the minimal 
5 = 2 state is disordered while all other states are ordered. On the frustrated pyrochlore lattice, 
we find (conservatively) that several states starting with the minimal S — 3 state are disordered. 
The disordered AKLT models we report here are a significant addition to the catalog of magnetic 
Hamiltonians in d = 3 with ground states known to lack order on account of strong quantum 
fluctuations. 



I. INTRODUCTION 

Quantum antiferromagnets have been a fertile field 
of research for a half century, exhibiting a great rich- 
ness and variety of physical phenomena. In more re- 
cent decades, starting with Anderson's introduction of 
the RVB stated and accelerating with the discovery of the 
cuprate superconductors^, much attention has focused on 
antiferromagnets that allow for disordered ground states 
due to a mix of frustration and quantum fluctuations'*. 
In an important step, Affleck et. alX showed how to con- 
struct models that build in a great deal of both these 
effects by using local projectors - models for which (essen- 
tially unique) ground states can be determined analyti- 
cally. These AKLT models have spins given \yy S — §Af , 
where M is any integer, and z the lattice coordination 
number. The associated ground states have the added 
feature that their wavefunctions can be written in Jas- 
trow (pair product) form. A general feature of such wave- 
functions is that the ground-state probability densities 
can be viewed as Boltzmann weights corresponding to a 
local, indeed nearest neighbor, Hamiltonian for classical 
spins on the same lattice. Using this unusual quantum- 
classical equivalence we can understand many properties 
of the states via Monte Carlo simulations of the associ- 
ated classical model. 

In d = 1 and d — 2, the AKLT states are disordered 
for any spin due to the Hohenberg-Mermin- Wagner the- 
orem. In particular, the d = 1 case is the celebrated 
AKLT chain which realizes the 5=1 Haldane phase. In 
this paper, we study AKLT states in d = 3, which are 
relatively less well-understood than their one and two di- 
mensional counterparts. Moreover, in three dimensions, 
the Hohenberg-Mermin- Wagner theorem no longer ap- 
plies and therefore whether an AKLT state of a given spin 



is disordered or instead exhibits long range order is no 
longer automatic. Instead a computation is now required 
to settle this question and it is this issue that we address 
in this paper by a combination of mean-field arguments 
and Monte Carlo simulation. Specifically, we discuss the 
AKLT states on the simple cubic and diamond lattices, 
where there is no (geometrical) frustration, as well on 
the highly frustrated pyrochlore lattice, where the atten- 
dant complications lead to a macroscopic ground state 
degeneracy of the associated classical model. Of course, 
all the models we study have frustration from competing 
interactions. 

On the cubic lattice we find that all AKLT states start- 
ing with the "minimal" (smallest spin) S = i state are 
ordered with the standard two sublattice Neel pattern. 
The diamond lattice has a small coordination number 
and thus larger fluctuations and we find that on it the 
the minimal S = 2 state is disordered while all higher spin 
states are ordered with the two sublattice Neel pattern. 
On the pyrochlore lattice, the geometrical frustration of 
the lattice plays a significant role. In mean field theory 
for the companion classical model we find a macroscopic 
number of solutions corresponding to as many energy 
minima. While the mean field estimate for the critical 
spin (transition temperature) already indicates that the 
minimal 5 = 3 model on the pyrochlore lattice is dis- 
ordered, the large number of competing states indicate 
that the true boundary between disorder and some form 
of order lies at much larger values of spin. Indeed, a 
basic simulation leads to a conservative bound in which 
disordered ground states persist up to S — 15. Given 
the unphysical complexity of the AKLT Hamiltonians at 
such large spins we do not pursue a more precise deter- 
mination of this boundary in this work. Indeed, readers 
may take as the main fruits of our work the identifica- 



tion of the S ^ 2 AKLT model on the diamond lattice 
and the S = 3 AKLT model on the pyrochlore lattice 
as (not too common) instances of three dimensional spin 
Hamiltonians with quantum disordered ground states. 

It is worth noting that models with quantum disor- 
dered ground states are currently objects of intense in- 
terest in the context of topological order and more specif- 
ically, in the context of topological quantum computing. 
We note that our disordered models do not yield topo- 
logically ordered states; they do not posess a topological 
degeneracy or host fractionalized excitations. The disor- 
dered states herein are described either as fully symmet- 
ric valence bond solids or, in the long wavelength sense, 
as quantum paramagnets. To understand why this is 
the case, it is instructive to recall how a closely related 
strategy works to produce topologically ordered states 
in S* = 1/2 models, including instances in d = 3. This 
strategy, initiated by Chayes, Chayes and Kivelson^, and 
brought to fruition in work by Raman, Moessner and 
Sondhi^ works with spin-i analogs of the AKLT mod- 
els called Klein models^. Unlike AKLT models, Klein 
models have many ground states — indeed they select 
the macroscopically many nearest neighbor valence bond 
coverings of a lattice. This selection of a degenerate man- 
ifold underlies the emergence of topological order. More 
precisely, the work of RMS showed that Klein models 
could be controUably perturbed on a family of lattices 
in order to select a topologically ordered (RVB) state in 
this ground state manifold. In this fashion they could 
construct SU(2) symmetric models with Z2 topological 
order in d — 2 but also models with Z2 and U(l) order 
in d = 3^. 

The rest of this paper is organized as follows: In Sec- 
tion II, we present a brief summary of the AKLT con- 
struction. We then proceed in Section III to review the 
mean field analysis of the AKLT states^. We then spe- 
cialize in Section IV to bipartite lattices, and compute the 
transition temperature for the simple cubic and diamond 
lattices using Monte Carlo simulations. We determine 
that while the simple cubic lattice exhibits Neel order for 
all choices of M (and thus S), the diamond lattice allows 
a quantum disordered state in the Af = 1 (5 = 2) case. 
We then go on to discuss the AKLT states on the frus- 
trated 3D pyrochlore lattice, and discuss the mean field 
analysis and classical ground states in this case. We find 
that the pyrochlore lattice admits quantum disordered 
states for many values of M; while the exact value of Mc 
was not determined, we find evidence from simulations 
that it exceeds 5, corresponding to S" = 15. 



II. AKLT STATES: A BRIEF REVIEW 

The central idea of the AKLT approacbi, is to use 
the idea of quantum singlets to construct correlated 
quantum-disordered wavefunctions, which are eigen- 
states of local projection operators. One can then pro- 
duce many-body Hamiltonians using projectors that ex- 



tinguish the state, thereby rendering the parent wave- 
function an exact ground state, typically with a gap to 
low- lying excitations. A general member of the family of 
valence bond solid (AKLT) states can be written com- 
pactly in terms of Schwingcr bosons"'^- : 
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This assigns M singlet creation operators to each link 
{ij) of a lattice C . The total boson occupancy per site is 
given by zM, where z is the lattice coordination number, 
and the resultant spin on each site is given by S" = ^zM. 
Thus, given any lattice, the above construction defines a 
family of AKLT states with S — ^zM, where M is any 
integer. The maximum possible spin on any link is then 
5,J°^'^ = 2S-M, and therefore |*(/:; M)) is extinguished 

by any Hamiltonian constructed out of projectors Pj{ij) 
onto link spin J, provided 2S - M + 1 < J <2S. The 
projectors, which transform as SU(2) singlets, may be 
written as polynomials in the Heisenberg coupling S'j • 5 ■ 
of order 25. Explicitly, one has 
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The AKLT states have a convenient representation in 
terms of SU(2) coherent states, as first shown in Ref . UXJ. 
In terms of the Schwinger bosons, the normalized spin- 
s' coherent state is given by \n) = (p!)^-^/^ (^^j^J,)^ I 0), 
where p = 2S, with z = {u , v) a CP spinor, with 
u = cos(6'/2) and v — sin(0/2)e"^. The unit vector 
n is given by n" = zV'^z, where a are the Pauli ma- 
trices. In the coherent state representation, the gen- 
eral AKLT state wavefunction is the pair product ^ = 
ri/ri) (^i ^j ~ "^i "j)*^- Following Ref. ^Q, we may write 
'I'P = exp(— i/ci/r) as the Boltzmann weight for a clas- 
sical 0(3) model with Hamiltonian 
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at temperature T = 1/M . All equal time quantum corre- 
lations in the state l^") may then be expressed as classical, 
finite temperature correlations of the Hamiltonian H^^. 

The consequences of this quantum-to-classical equiva- 
lence, which is a general feature of Jastrow (pair prod- 
uct) wavefunctions, were noted in Ref. 'lO'. On one 
and two-dimensional lattices, the Hohenberg-Mermin- 
Wagner theorem precludes long-ranged order at any finite 
value of the discrete quantum parameter M . Thus, while 
the S — 2 Heisenberg model on the square lattice is rig- 
orously known to have a Neel ordered ground stateii, the 
S — 2 AKLT Hamiltonian, which includes up to biquar- 
tic terms, has a featureless quantum disordered ground 
state, called a 'quantum paramagnet'. In three dimen- 
sions, we expect Neel order for large M, corresponding 



to low temperatures in the classical model. If the Neel 
temperature for 7?^j on a given lattice satisfies T^ > 1, 
then Mj. < 1, and all the allowed AKLT states on that 
lattice exhibit long-ranged order. 

The issue of whether or not the AKLT states can be in 
the quantum disordered phase on a given lattice can be 
investigated by a combination of mean-field calculations 
and classical Monte Carlo simulations, which we present 
below. 



III. MEAN FIELD THEORY 



ignores fluctuations, hence it overestimates T^ and un- 
derestimates Mj.. Therefore the possibility remains that 
a quantum disordered AKLT state may exist in a three- 
dimensional lattice. We examine two cases, the simple 
cubic lattice (z = 6) and the diamond lattice (z = 4). 
We shall address this issue via classical Monte Carlo sim- 
ulations of our model on both lattices. We note that, of 
these, the diamond lattice is the stronger candidate as it 
is more weakly coordinated, and M"'' = | is sufficiently 
close to threshold that fluctuations are likely to drive the 
true M^ to be greater than unity. 



A mean field analysis of the classical model of eqn. 
13] on bipartite lattices was described in Refs. MUX In 
the general case, we may begin with the Hamiltonian of 
eqn. [31 and we write n^ = m^ + Sh^, with (n^) — rfi^. 
Expanding iJ^, to order 5n^, we obtain the mean field 
Hamiltonian H^^^ = Eq — Tlii^i ' ^ii where the mean 
field /ij is given by 



— rrii ■ mi 



(4) 



where the prime restricts the sum on j to nearest neigh- 
bors of site i. Self-consistency then requires 

J '^n^ exp (/i, • nJT) 
which yields fh^ = m^ h^j\h^^ with local magnetization 



m. = coth I ^) - ^ • (6) 



IV. UNFRUSTRATED LATTICES: SIMPLE 
CUBIC AND DIAMOND 

A. Mean Field Transition 

On an unfrustrated, bipartite, lattice a sublattice ro- 
tation fij ^ rj^h^, with r/j — ±1 on the A (B) sublattice, 
results in a ferromagnetic interaction, and if we posit a 
uniform local magnetization m we obtain the mean field 
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where z is the lattice coordination number. This results 
in a mean field transition temperature T^^ = ^z, i.e. 
M^" = 3^"^ AU AKLT states on bipartite lattices in 
more than two space dimensions will exhibit two sub- 
lattice Neel order, provided M > M^. According to 
the mean field analysis, long ranged order should per- 
tain for z > 3, which would be satisfied by almost any 
three-dimensional structure. However, mean field theory 



B. Monte-Carlo Simulations 

The classical Hamiltonian H^^ of eqn. [3] consists 
of nearest neighbor interactions v('^^A where t9j = 
cos^^ (rij-n) is the relative angle between spins on neigh- 
boring sites i and j, and v{i9) — — lnsin^(^i?). This 
interaction strongly suppresses ferromagnetic alignment, 
with a logarithmically infinite barrier, but has a smooth 
quadratic minimum v{'d) w 1(19 — tt)^ when d k, tt. We 
have simulated the equivalent ferromagnetic model, with 
interaction v{d) — — lncos^(|^). 

We used a multithrcad Monte Carlo approach, in which 
simultaneous simulations with independent initial config- 
urations were used to produce M independent Markov 
chains each with N configurationsi^, which were then 
written to a file. For every independent thread, we per- 
formed checkerboard sweeps of the lattice using a stan- 
dard Metropolis Monte Carlo technique^ -. In each Monte 
Carlo step, we produced a vector 5fi, with length dis- 
tributed according to a Gaussian and pointing in a ran- 
dom direction, which was used to generate a new spin 
unit vector 



5fi 



(8) 



The standard deviation of the Gaussian was adjusted by 
hand until a significant fraction of proposed moves were 
accepted (we left it at a = 0.5.) 

The number of Monte Carlo steps per site (MCS) 
and the number of independent threads were adjusted to 
count roughly the same number of autocorrelation times 
for each sample size^^. For each chain, we obtained the 
average value of the Binder cumulant, and averaged this 
across chains to get a single number for each tempera- 
ture. We estimated the error from the standard devia- 
tion of the M independent thread averages. This is free 
of the usual complications of correlated samples inherent 
in estimating the error from a single chain, and frees us 
of the need to compute autocorrelation times to weight 
our error estimate. 

Plots were made of the Binder cumulanti^, defined to 
be 
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where M = ^^ n^ is the total magnetization. 

For any system of Heisenberg spins in the thermody- 
namic Hmit, the Binder cumulant has value | in the low- 
temperature (ordered) phase and value g in the high- 
temperature phase. These are easily seen by assuming a 
gaussian distribution for \M\ at high temperature, and 
using the result that all the expectation values of powers 
of M ■ M are equal in the ordered phase. For a finite 
system, the limiting values continue to be close to these 
estimates, but the interpolating behavior is different for 
each system size; the primary utility from our point of 
view is that finite-size scaling analysis of B reveals that 
it has a fixed point at the transition tcmperatureiS. We 
may therefore determine Tc by plotting the Binder cumu- 
lant for a series of different lattice sizes, and determining 
the points where the curves cross. 

Before simulating our modified interaction, we checked 
our code by determining the (known) transition temper- 
atures for the standard Heisenberg model on the simple 
cubic latticei^ and the Ising model on both the diamond 
and the simple cubic lattices^^, as well as comparing 
the high-temperature susceptibility from simulations to 
the predictions of the high-temperature expansionii. All 
these agreed well with the expected values, at least to the 
accuracy we need to determine whether T^. is less than 
or greater than 1. Recall that if Tc < 1, then Mc > 1, 
which means that the minimal AKLT state, with M = 1, 
is on the disordered side of the phase transition. 

Using our Monte Carlo simulations, we obtain esti- 
mates of Tc for the families of AKLT states on the sim- 
ple cubic and diamond lattices. Although our simulation 
techniques were not particularly sophisticated, they were 
sufficient to pin down Tc to a reasonable degree of accu- 
racy, and certainly enough to determine whether Tc > 1. 
Our simulations allow us to estimate that Tc*^ « 1.66 
on the simple cubic lattice, and that T^ w 0.85 for the 
diamond structure. Therefore, we conclude that while 
all the simple cubic AKLT states are Neel-ordered, the 
minimal {S = 2) AKLT state in diamond is a featureless 
quantum disordered state. 



FRUSTRATED LATTICE: 
PYROCHLORE 



THE 



The pyrochlore is a lattice of corner-sharing tetrahe- 
dra and can be constructed from the the diamond lattice 
by placing a site at the midpoint of each bond, result- 
ing in a quadripartite structure. The pyrochlore lattice 
is highly frustrated from the perspective of of coUinear 
antiferromagnetism; the canonical nearest-neighbor clas- 
sical Heisenberg antiferromagnet on this lattice has an 
extensive ground-state degeneracy and remains a quan- 
tum paramagnet at all temperaturesi^. 

Our problem has a different form for the interaction 
and hence the results for the nearest neighbor problem, 
which build on the high degree of degeneracy for a single 
tetrahedron, do not apply. Indeed, as we discuss below. 
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FIG. 1: Binder cumulant plots for the valence-bond states on 
the cubic and diamond lattices. The T-axis scale is chosen 
using a rough estimate of Tc so as to provide approximately 
the same window in natural units T/T^ for both cases. In 
each case, the total number of spins being simulated is 2 ■ L^. 
We perform a fit of the data (weighted by the error bars) 
to a parabola and estimate Tc from the intersection of the 
best-fit lines. We can be reasonably confident that the curves 
have an intersection from the fact that they change order 
on either side of the crossing, and become separated by more 
than a standard deviation as we move away from the crossing. 
We obtain Tc ~ 1.66 for the cubic lattice and Tc ~ 0.85 for 
diamond. 



the logarithmic form of the interaction energy leads to the 
selection of a unique single-tetrahedron ground state up 
to global rotations. However, the full lattice still exhibits 
a substantial ground state degeneracy on account of its 
open architecture indicating an anomalously low transi- 
tion temperature which we roughly bound from above by 
T « 0.2. We now turn to the details of these assertions. 



A. Single- Tetrahedron Ground States 

Numerical minimization on a single tetrahedron finds 
the lowest-energy configuration to be one where each pair 




FIG. 2: The quadripartite pyrochlore lattice, which is formed 
out of corner-sharing tetrahedra. 



of spins make an angle dij = cos~^{h^-nj) — cos^^ (~i)- 
This means that the spins are pointing either towards or 
away from the corners of a regular tetrahedron in three- 
dimensional spin space. We proceed to search for soft 
modes, by expanding the energy to quadratic order and 
studying the resulting normal modes. We find that there 
is a pair of soft modes corresponding to global rotations, 
and another where three spins rotate about the axis de- 
fined by the fourth. The latter mode leads to a degener- 
acy of ground states of the full lattice as discussed below. 
We note that for the Heisenberg antiferromagnet with 
interaction n^ • fij , the single tetrahedron Hamiltonian is 

Hr = (^fr) I where Mr — J2ier "-« i^ ^ ^^™ ^^ ^^'^ ^Pi^^ 
vectors over all sites in the tetrahedron F. The ground 

state manifold Mr = is then five-dimensional, since 
one can choose any two vectors fiA and fiB, then take 
fic = —fiji and n]j = —fis- The four freedoms associ- 
ated with choosing fiA and fiB are then augmented by 
another freedom to rotate the C and D spins about the 
direction fiA + ns- A large- A^ analysis^^ finds that the 
0{N) pyrochlore antiferromagnet is paramagnetic down 
to T = 0. 



state Potts model. We therefore conclude that one family 
of ground states of the classical Hamiltonian on the py- 
rochlore lattice are in a one-to-one correspondence with 
the ground states of the 4-state Potts antiferromagnet on 
the pyrochlore lattice. Readers familiar with the lore on 
the kagome problem2£ will recognize the resemblance to 
the planar ground states there which are in correspon- 
dence to ground states of the 3-state Potts model. As in 
the kagome problem, from this set of ground states oth- 
ers can be constructed by identifying sets of spins which 
can be locally rotated by an arbitrary amount at zero 
energy cost. These are sets of spins, say of type B, C 
and D which are connected to other spins solely by spins 
of type A. Clearly one can rotate this set by an angle 
about the A axis at zero energy cost. 

While we have not parametrized the full, continuous, 
space of ground states an extensive lower bound on the 
degeneracy of the 'Potts submanifold' of ground states 
can be obtained as follows. First, we note that the 
number of allowed configurations of the 3-state Potts 
model on a kagome lattice with M sites is given^Sii^ by 
.gk ~ (1. 20872) (2AV3)^ ^^^^ ^g partition the pyrochlore 
into four sublattices, so that the sites that lie on a single 
tetrahedron are each on different sublattices. Choose one 
sublattice, and fix the spins on that sublattice to be one 
of the four types (say A.) Now, looking down through the 
tetrahedra, one sees alternating layers of triangular and 
kagome planes; the kagome planes are made up of B, C, 
and D spins, while the triangular planes are made up of 
A spins. In each kagome plane, we have M spins, whose 
configurations are those of the 3-state Potts model. If 
we now let A^k be the number of kagome planes, we must 
have that M ■ iVk = |A^, where N is the total number 
of spins in the system. We then have for the number of 
states in this restricted submanifold 
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(10) 



where the factor of 4 stems from the fact that we can 
choose any one of the four spins to be fixed in the trian- 
gular planes. Since we've restricted ourselves to consid- 
ering a certain submanifold of the ground states in the 
above argument, it is clear that we have obtained a lower 
bound for the degeneracy of the Potts submanifold. 



B. Ground States on the Full Lattice 

There are many ways in which we can construct degen- 
erate states on the lattice that simultaneously satisfy the 
minimum-energy constraint on every tetrahedron. We 
begin by describing the simplest such states which form 
a discrete family. To this end, label the four spins de- 
fined by the single-tetrahedron constraint (with a fixed 
joint orientation) as A, B, C and D. If we use only 
these four orientations for each spin, we have the con- 
straint that none of them can occur twice on the same 
tetrahedron; this translates to the statement that spins 
on neighboring links must be different. This is the same 
constraint as for ground states of the antiferromagnetic 4- 



C. Bounds on Tc 

Each of the ground states identified above can serve as 
a basis for a mean-field treatment of the system and all of 
them yield the same T^^. This vast set of "soft modes" 
is, of course, a signature that the true Tc <C T"''. Thus 
we may begin with a calculation of T"'' which can serve 
as an upper bound on the true Tc- 

Consider a spin at site i in the pyrochlore lattice. Ex- 
panding in small fluctuations about any ground state, 
we have the same mean-field Hamiltonian as in the gen- 
eral mean field Ansatz of section IIIIl with the mean 
field given by eqn. 3] In a mean-field treatment, each 



of the neighbor spins fij is to be replaced by its aver- 
age ruj = {fij) — mej in the particular ground state 
that we are considering. In any ground state, we note 
that the angle between any pair of nearest neighbors is 
i^ij = cos^^ (~i)- I^ addition, the spins on a tetrahe- 
dron add to zero, which allows us to write ^- rfij = —mi. 
If we further recall that each spin lies on exactly two 
tetrahedra, we obtain the following expression for the 
mean field acting at site i: 



K = - 
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(11) 



Note that we have only made use of the local structure 
of the ground state, and so our treatment here is rele- 
vant for the transition into any state in the ground state 
manifold. Substituting the mean field in eq. [11] into the 
self-consistency condition (eqnl6]), we find, in a manner 
similar to the bipartite case, that the mean-field esti- 
mate of the transition temperature is T"'' = -^. From 
this alone we conclude that the M = 1 state on the py- 
rochlore lattice is quantum disordered. 

There is little question that the actual T^ is much 
lower than the mean-field estimate and therefore Mc is 
much higher than | , allowing many more quantum disor- 
dered states. As is familiar from other highly frustrated 
magnets, where the ground-state manifold encompasses 
a vastly degenerate set of states the transition will be 
driven by the 'order-by-disorder' mechanism wherein a 
particular state or subset of states is favored by entropic 
effects at low temperatures. This is a weak effect and 
hence Tc is typically a small fraction of T^^^ . From coarse 
Monte Carlo simulations we find evidence that Tc < 0.2 
corresponding to Sc > \zM — 15. This leads already 
to an AKLT Hamiltonian that is a degree-60 polynomial 
in the spins and thus there is little reason in the cur- 
rent context to locate the transition or the nature of the 
ordered phase with greater precision. 

However, we note that the same set of ground states 
arises in the the classical Heisenberg model with nearest 
neighbor bilinear and biquadratic interactions with the 
latter chosen to disfavor collinearity. This is a physically 
plausible model and we will report a fuller investigation 
of it elsewhere^^. 



VI. CONCLUDING REMARKS 

To summarize, we have studied AKLT states on two 
unfrustrated and one frustrated lattice in d = 3 by a 
combination of mean-field theory and Monte Carlo sim- 
ulations for the associated classical models. Wc find that 
the simple cubic lattice is Neel ordered at all values of 
the singlet parameter AI and spin S; the diamond lat- 
tice, on the other hand, is quantum-disordered for AI — 1 
(5 = 2), and Neel ordered for Af > 1. On the pyrochlore 
lattice we find that the M = 1 (5 = 3) model is definitely 
disordered and the boundary between disorder and order 
very likely lies above M = b. 

While quantum-disordered ground states in low (i.e. 
one and two dimensions) have often been discussed, three 
dimensions has historically been the province of long 
range order. Hence our disordered models on the di- 
amond and pyrochlore lattices significantly expand the 
set of possibilities for quantum ground states of models 
with Heisenberg symmetry in d = 3. 

In recent work, one of us has generalized the AKLT 
construction to SU(A^) spins^. In the near future, we 
intend to investigate the SU(4) simplex state on the py- 
rochlore lattice introduced in this work by methods sim- 
ilar to the ones used in the present paper—. 
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